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Abstract 

We study properties of the lightest neutralino (%°) and calculate its cosmological relic density in 
a supersymmetric U(l)' model with a secluded U(l)' breaking sector (the /S-model). The lightest 
neutralino mass is smaller than in the minimal supersymmetric standard model; for instance, 
m x o ^ 100 GeV in the limit that the U(l)' gaugino mass is large compared to the electroweak 
scale. We find that the Zx°X° coupling can be enhanced due to the singlino components in the 
extended neutralino sector. Neutralino annihilation through the Z-resonance then reproduces the 
measured cold dark matter density over broad regions of the model parameter space. 
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I. INTRODUCTION 



Recent mapping of the Cosmic Microwave Background (CMB) anisotropy [l| has provided 
precision information on the densities of matter and dark energy in the Universe. The major 
part of the matter is non-relativistic and non-baryonic (cold and dark). When the Sloan 
Digital Sky Survey (SDSS) data on large scale structure are analyzed in combination with 
the Wilkinson Microwave Anisotropy Probe (WMAP) data, a cold dark matter (CDM) relic 
density 

n CDM h 2 = 0.12 ±0.01 (SDSS + WMAP) (1) 

r—i 

is found [2j. This very restrictive range for flcDMh has significant impact on the allowed 
masses and couplings of cold dark matter particles. 

A neutral, stable, massive particle that interacts weakly is a natural candidate for CDM. 
The Standard Model (SM) does not have a particle with these properties but a supersym- 
metric model with i?-parity conservation, one of the best motivated new physics possibilities, 
does. In the Minimal Supersymmetric Standard Model (MSSM) the lightest neutralino is 
the favored lightest supersymmetric particle (LSP). 

The lightest neutralinos existed abundantly in the early Universe in thermal equilibrium 
with other particles, where their pair annihilations were balanced by pair creation. As 
the Universe cooled, the neutralino density became Boltzmann-suppressed. Deviation from 
thermal equilibrium began when the temperature reached the freeze-out temperature Tf ~ 
m x o /20. After the temperature dropped to ~ |Xy, the annihilation rate became equal to the 
expansion rate, and the neutralino relic density was n x o = H/ (crv), where H is the Hubble 
expansion rate at that temperature. Here (crv) is the thermally averaged cross section times 
neutralino velocity 3]. The remaining CDM density relative to the critical density is 

Q x o = n x om x o/p c = Hm x o/ ((crv) p c ) . (2) 

A viable CDM model must reproduce the recent precise measurement of I^cdm, preferably 
without fine-tuning of the model parameters. The full test of neutralino dark matter can 
be accomplished with the direct detection of the LSP in collider experiments and/or in 
elastic scattering experiments at underground detectors. The parameters of the model must 
quantitatively explain Q x oh 2 ~ 0.12 and the direct detection rates. 

There have been extensive theoretical studies of the relic density 0, S 0] • In the minimal 
supergravity model (mSUGRA) annihilations of stable supersymmetric particles can repro- 
duce the right order of magnitude. However, the allowed regions of mSUGRA parameter 
space are becoming increasingly constrained jjj by the combination of the recent relic den- 
sity and the LEP data and now allow only a few regions of parameter space including (i) 
pair annihilation of neutralinos with dominant bino composition through A , if - resonances 
at high tan/5 (the so-called A- funnel) jH, (ii) neutralino-stau coannihilation and (iii) 
annihilation of neutralinos with mixed gaugino-Higgsino components in the focus point re- 
gion 0| . A nearly pure bino LSP state can give the right size of the relic density 
and also satisfy the required radiative electroweak symmetry breaking (EWSB) in unified 
models 13. In SO (10) grand unified supergravity models with Yukawa unification, it is 
more difficult to obtain the required relic density, but it is still possible with specific scalar 
mass patterns 
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The MSSM suffers from the /^-problem |14J and the lack of a sufficiently strong first order 
phase transition for electroweak baryogenesis (EWBG) over most of the parameter space 1 . 
The missing ingredients of possible TeV-scale new physics that would cure these problems 
may modify the properties of the CDM candidate. The Next-to-Minimal Supersymmet- 
ric extension of the Standard Model (NMSSM) [16] avoids the problems by introducing 
an additional SM singlet Higgs and a discrete symmetry Z 3 , but it then suffers from the 
cosmological domain wall problem [l7| . 

A natural extension of the MSSM that avoids the above difficulties is a supersymmetric 
U(l)' model, with Higgs singlets to break the additional Abelian symmetry spontaneously 
at the TeV scale 2 . Additional U(l) gauge symmetries are predicted in many types of new 
theories, including superstrings [20|, grand unification, extra dimensions 21], dynamical 
symmetry breaking [22j, and the little Higgs model [23]. In this letter we investigate the 
properties of the lightest neutralino and evaluate its relic density 3 in an extended model of 
the MSSM with an extra U(l) gauge symmetry and four extra Higgs singlets (S, Si, S 2 , and 
S3) [H|. The superpotential is 4 

W = h a SH x H 2 + X s SiS 2 S 3j (3) 

where h s and A s are dimensionless parameters. We call this the 5-model. 
The .S-model solves the /i-problem by generating an effective fi parameter 

fi cS = h s (S) . (4) 



The model is free from the domain wall problem since there is no discrete symmetry |2j 
I2R I27I ] . The .S-model can also provide a sufficiently strong first order phase transition and 
new sources of CP violation for EWBG [2^ . The Z' has a large mass 5 generated by singlet 
Higgs fields, 61,2,3, which acquire large (TeV scale) VEVs for small X s because of an almost 
F and D flat direction. These multiple singlets allow /z e ff to be of the electroweak scale while 
keeping the Z' heavier than the experimental limit. 

Electroweak symmetry breaking is driven by electroweak scale trilinear terms. This leads 
to tan/3 = v%/vx ~ 1, while solutions without unwanted global minima at (Hf) = typically 
have (S) ^1.5 (Hf) jiij]. However, both of these conditions can be relaxed somewhat |31j . 
All dimensionful supersymmetry-breaking parameters are close to the electroweak scale. The 
squark and slepton masses are similar to those of the MSSM, and the soft supersymmetry 
breaking Higgs and singlet masses are smaller than the Higgs doublet mass in the MSSM. 

In our analysis, we consider the limit with M[ 3> Mi^, where M[, M\, and M 2 are the 



gaugino masses for U(l)', U(l)y, and SU(2)l, respectively [2Jj|3l|]. This limit suppresses 



the effects of the U(l)' Higgs charges and greatly reduces the number of free parameters 



1 To have a strong first-order phase transition in the MSSM, the light Higgs mass should be only slightly 
above the LEP experimental bound and the light stop should be lighter than the top fl5| . 

2 For another extension of the MSSM using a discrete symmetry but free of domain wall problem, see the 

nMSSM model (HQ. 

3 For an earlier study of the relic density in a U(l)' model with one singlet, in a different approach and 
framework, see Ref. [2^. 

4 For a U(iy model with only one singlet, see Ref. |26(. It is more difficult to obtain Mz> ^> Mz in such 
models. 

5 The CDF limit is M z > ^ 500 - 800 GeV, depending on the model 0,13(1 . 
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in the neutralino mass matrix. The allowed mass range of the lightest neutralino then is 
limited (m x o < 100 GeV) compared to that of the MSSM, due to the large component of 
the singlino in the neutralino, which we assume is the LSP. While the small mass of the LSP 
makes the heavy fermion and gauge boson channels of the MSSM irrelevant, it opens up new 
relevant channels, particularly the Z-resonance and the light Higgs resonance channels. The 
lighter Higgs scalar (h°) in this model is a mixture of the Higgs doublets and singlets, and 
it can be much li ghte r than the LEP bound of m h o > 115 GeV that applies to the SM-like 
Higgs bosons |25l.l3ll] . 

In the MSSM, Z-resonance annihilation is not likely to be a relevant channel because the 
Zx°X° coupling is small when tan/5 ~ 1. However, when the singlino component is intro- 
duced into the neutralino sector, the Zx°X° coupling can be enhanced to give sufficiently 
large x°X° annihilation. This Z-resonance annihilation alone, with suitable parameter val- 
ues, can reproduce the acceptable cold dark matter relic density in most of the allowed x° 
mass range in the S'-model. 



II. MASS AND COUPLING OF THE LIGHTEST NEUTRALINO 

We consider a scenario with a massive Z' and VEVs of Si (i = 1,2,3) that are large 
compared to other electroweak scale parameters. There is an approximate decoupling of the 
neutralinos associated with the Z' and the Si, and the effective neutralino mass matrix for 
the remaining neutralinos in the basis of j-B, W3, H®, H®, S^ is 6 



( Mi 

M 2 

-giVx/2 g 2 vi/2 

giv 2 /2 -g 2 v 2 /2 





-gm/2 
g 2 vi/2 


-h s s/V2 
-h s v 2 /V2 



g x v 2 /2 

-g 2 v 2 /2 

-h s s/V2 


-h s Vx/y/2 






-h s v 2 /V2 
-h s vxl\[2 




(5) 



where e = g± cos 0w = g 2 sin 9 



w ■ 



% with 



The VEVs of the Higgs doublets are (H®) 

246 GeV and the VEV of the Higgs singlet is (S°) = -j=. The mass eigenstates 
are ordered as m y o = m x o < m x o < ■ ■ ■ < m x o . This mass matrix leads to a kind of see-saw 
mechanism (lil EH l33| that makes the lightest neutralino mass very small. In our analysis 
we scanned over the parameters 



tan/3 = 1.03,1.5,2.0,2.5 

h s = 0.1 to 0.75 (in steps of 0.01) 
M 2 = -500 GeV to 500 GeV (in steps of 1 GeV) 
s = 100 GeV to 500 GeV (in steps of 1 GeV) 

p. 2 

We use the gaugino mass unification relations (Mi = |%M 2 ~ 0.5M 2 ), but take M[ > M 2 . 
The U(l)' charge dependence vanishes from the effective mass matrix for large M{. 



6 In this limit, the neutralino matrix is basically the same as that of the one-singlet models and shares 
many of the same properties Q, ■ 



4 



The MSSM mass matrix 7 corresponds to dropping the last row and column in the matrix 
of Eq. © and taking \i = h s s/V2. The MSSM upper bound on m x o is very sensitive to 
the value of M 2 , while that of the S-model depends sensitively on the value of h s . In both 
models, the upper bound on the lightest neutralino mass has its largest value when tan/? 
is 1.03 (or l) 8 . When tan/5 ~ 1, the m x o bound in the S'-model increases almost linearly 
with h s . We choose our upper bound of h. to be 0.75 as advocated in Ref. (Hf to keep the 
running of h s finite up to Plank scale 0,0. 

We find that 

m x o < 99 GeV. (6) 

The maximum value of m x o occurs with tan/3 = 1.03, h s = 0.75, M 2 = —183 GeV and 
s = 130 GeV for the S'-model 9 while for the MSSM the maximum m x o < 254 GeV occurs 
for tan/5 = 1.03, M 2 = -500 GeV and /i ~ 261 GeV. 

The small LSP mass makes the Z- resonance contribution to the relic density relevant in 
the S-model. For Zx^Xj and Z'xiX^j interactions, the general Lagrangian is given by 

c = -a z z, (Hh^m - mi^m 



+ (2Q 1 if 1 V 7 5 #i + 2Q 2 if 2 V7 5 #2° + 2Q 5 SV 7 5 S) (7) 

which, in terms of neutralino mass eigenstates, can be written as 

tzvw = \9zZ^ (ofP L + 0'fP R ) x ° 

+ 19z'Z^Y(i^Pl + I^Pr)xI (8) 
Here the couplings are defined as 

°lf = ~I n *n* 3 + \n u n; a of = ±n; 3 n 33 - l -N*,N 3i (9) 

if = -QxN i3 N* 3 - Q 2 N li N* A - Q s N i5 N* 5 (10) 
if = Q 1 N^N j3 + Q 2 N* i4 N J , + Q s N* 5 N j5 , (11) 

where the matrix N relates the weak and mass bases. Qx, Q 2 and Qs are the U(l)' charge 
of the H%, # 2 ° and S, respectively. The Z X °X° coupling 0'[[ = -0'{? = -\ (|A^ 13 | 2 - |7V 14 | 2 ) 
is composed of the two doublet Higgsino components coupling to the Z boson. A massive 



7 The MSSM parameter range considered here (especially for tan /3) is not experimentally allowed (e.g. by 
the lightest Higgs mass), but is chosen to demonstrate the effects of adding the U(l)'. 

8 Since the exact value tan/3 = 1 would make the coupling of Zx°X° vanish in both the MSSM and the 
S'-model, we choose tan/3 = 1.03 instead of 1. 

9 The upper bound on m x o in the different limits are: (i) for M[ ^ M 2 , Si ~ O(EW), m x o < 100 GeV; 
(ii) for M{ = Mi, s z > O(EW), m x « < 280 GeV; (iii) for M{ = M u s z ~ O(EW), m x o < 100 GeV. We 
scanned 100 to 1000 GeV for Mi and Si and assumed h s < 0.75 with A s < 0.2. For the U(l)' couplings, 
various Eq model charge assignments (x^^v) an( i the GUT motivated coupling constant gz* = yfs'i 
were assumed. For a full 9x9 neutralino mass matrix, see Ref. 25l Kill. 
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FIG. 1: The fraction of the lightest neutralino (|iVij| 2 ) versus M 2 for (a) the S'-model and (b) 
the MSSM. All components are presented: \Nu\ 2 (B, dot), |iV 12 | 2 (W3, dash), | iV 13 1 2 (Hf, solid), 
|iVi4| 2 (H°, dash-dot), and \Ni 5 \ 2 (S, dash-dot-dot). The curve for |Aii| 2 is labeled by i. Fixed 
values of h s = 0.75, tan/3 = 1.03 and s = 250 GeV are illustrated. The difference of | A^i3 1 2 and 
\Ni4\ 2 is large when the singlino component is present even though tan/5 ~ 1. 

(TeV-scale) Z' would not provide sufficiently large annihilation for our light neutralinos, and 
thus in our analysis we use only the Z-resonance annihilation. 

FigureJU presents the 1 7V Xi | 2 with i = 1 (B, dot), 2 (W 3 , dash), 3 (H$, solid), 4 (H$, dash- 
dot), 5 (S, dash-dot-dot), for M 2 = —500 to 500 GeV and a fixed set of values h s = 0.75, 
tan/3 = 1.03 and s = 250 GeV. The value of each |iVij| 2 is shown versus M 2 in (a) the 
S-model and (b) the MSSM. We note that in the region around M2 ~ 300 GeV the singlino 
component |iV 15 | 2 increases, and |iV 13 | 2 deviates substantially from (A^ 2 making the Zx°X° 
coupling large. This does not happen in the MSSM for the same parameters. In addition, 
for M2 —150 GeV, the singlino dominates and the difference of | iV"i3 1 2 and | iVi^ | 2 in the 
S'-model is larger than that in the MSSM. As an example, the components of the LSP for 
M 2 = 330 GeV are 

X° = 0.185 - O.IIW3 - 0.09^° - 0.58# 2 ° + 0.78S (S-model), (12) 
X° = 0.465 - 0.27iy 3 + 0.60iJ 2 - 0.60^3° (MSSM). (13) 

The large difference in | iV"i3 1 2 and | A^i^ | 2 for the S'-model is remarkable in view of the fact 
that tan (3 = 1 makes the H® and if 2 parts in the mass matrix of Eq. (J3J) the same up to the 
sign, which leads one to expect |Ai 3 | 2 — |Ai 4 | 2 ~ 0, as is the case in the MSSM (Figure d(b)). 
Thus, the addition of the singlino component plays a critical role in enhancing the difference 
of |A^ 13 | 2 and |A^4| 2 , which is particularly important for the x° m t ne S-model to generate 
sufficient annihilation mediated by the Z boson. The suppressed coupling (by tan/5 ~ 1) 
can still be large enough to give an acceptable relic density. Figure Q also shows that the 
gaugino components |A^n| 2 (bino) and |A^ 12 | 2 (wino), especially the bino, are dominant for 
relatively small values of |M 2 | in both models. 

We present {O'^l = \0'{£\ = \0'{?\ for the S-model and the MSSM in Figure H For the 
same values of h s = 0.75 and s = 250 GeV, we select 3 different tan/5 values of (a) tan/5 = 
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FIG. 2: The Z X °X° coupling \0'{ x \ versus M 2 for (a) tan/5 = 1.03, (b) tan/3 = 2.0, and (c) 
tan/3 = 2.5, in the S-model (solid) and in the MSSM (dash). Fixed values of h s = 0.75 and 
s = 250 GeV are used. The coupling is much larger in the S-model than in the MSSM in most of 
the parameter space. 



1.03, (b) tan/3 = 2.0, and (c) tan/3 = 2.5. For each tan/?, the Zx°X° coupling is much larger 
in the S'-model than in the MSSM in most of the parameter space. We numerically checked 
that this feature holds for other choices of the parameter values (0.4 < h s < 1.0, 100 GeV 
< s < 1000 GeV, 0.5 < tan/3 < 2.5). 

For tan/3 ~ 1 (Figure 121 (a)), the Zx°X° coupling in the 5-model is at its maximum for 
M 2 ~ 300 GeV and there is a relatively small but still noticeable peak in the MSSM at 
M 2 — —250 GeV. The x° rei ic density depends on not only the coupling but also whether 
m x o is near Mz/2. 




FIG. 3: (a) The lightest neutralino mass as a function of M 2 and (b) the Zx X ° coupling (|Oii|) 
versus the lightest neutralino mass in the S'-model (solid) and the MSSM (dash). Fixed values of 
h s = 0.75, tan/3 = 1.03 and s = 250 GeV are used. The 5-model has a smaller m x o bound, and, 
for M 2 > 0, larger Z X ° X ° coupling than the MSSM. 
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Figure El shows (a) the x° mass versus M 2 and (b) \0'{ x \ versus the LSP mass. The same 
parameter values as Figures ^ and El (a) are used. The LSP mass in the S-model is almost 
constant for M 2 < -200 GeV and M 2 > 300 GeV, and it is smaller than in the MSSM. 
The lightest neutralino masses in both models have practically identical dependence on M 2 
before they reach the flat curves in our parameter choice (tan/3 ~ 1). Enhancements of 
couplings in the S'-model and the MSSM are found around m x o ~ 80 GeV and 130 GeV, 
respectively, that is, for the flat parts of the mass curves in Figure |3] (a). 

III. NEUTRALINO ANNIHILATION MEDIATED BY THE Z BOSON 

The relic density of the lightest neutralino is inversely proportional to the annihilation 
cross section of x° pairs. When kinematically allowed, the neutralino pairs annihilate into 
pairs of fermions, gauge bosons, and Higgs bosons through s, t, and u channel diagrams. 
In general, the annihilation cross section is greatly enhanced by the Z boson or the Higgs 
boson poles, and the relic density is correspondingly suppressed when 2m x o ~ Mz or 
2m x o ~ M,/,. 

The Higgs masses and couplings in the S'-model depend on several free parameters. More- 
over, the light Higgs width is very narrow and a fine-tuning of parameters would be required 
for neutralino annihilation through 0° to generate the appropriate relic density, and the 
other Higgs bosons would often be too heavy to have resonance effects because the LSP is 
light [3jJ. We are therefore interested in the case that the Zx°X° coupling is sizable and 
Z mediated annihilation alone can lead to an acceptable relic density (£7cdm^ 2 )- A more 
complete investigation is in progress to include all Higgs resonance effects. We require a 
lower bound on m x o of Mz/2 so that Z — > x°X° does not affect the Z width significantly 10 . 
Then the lightest neutralino mass in the S-model is constrained to the range 

46 GeV < m x o < 100 GeV. (14) 

with the upper bound from Eq. (0). 

We evaluate the annihilation cross section including interferences among all diagrams 
by calculating the amplitude of each diagram with the helicity amplitude formalism, then 
numerically evaluate the full matrix element squared. The neutralino relic density is cal- 
culated with relativistic Boltzmann averaging, neutralino annihilation threshold effects and 
Breit-Wigner poles (H El IH . 

We show the neutralino relic density (Q x oh 2 ) versus m x o in Figure 0] for both the S'- 
model 11 and the MSSM for the tan/3 choices (a) tan/3 = 1.03, (b) tan/3 = 2.0, and (c) 
tan/3 = 2.5. We included all possible channels for the MSSM calculation; for the S-model 
we included only the kinematically relevant s-channel processes with the Zx°X° coupling, 
X°X° -> Z -> fifi, W + W~ (fi is a SM fermion). We scanned over -500 GeV < M 2 < 500 
GeV for h s = 0.75, and s = 250 GeV to calculate m x o and evaluate the Zx°X° coupling. 
There are several interesting aspects to note. 



10 Since AFi nv < 2.0 MeV (95% C.L.) I3J], smaller m x o is possible [l9|. For reasonably smaller masses, our 
results are unchanged. 

11 For similar relic density results in the nMSSM, see Ref. [l9|. For the NMSSM relic density calculations, 
see Ref. 0. 
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FIG. 4: The neutralino relic density (Q x oh 2 ) versus m x o in the .S'-model [M2 > (solid) and 
M 2 < (dash-dot)] and the MSSM [M 2 > (dash) and M 2 < (dot)] for (a) tan/3 = 1.03, (b) 
tan/3 = 2.0, and (c) tan/3 = 2.5. The isolated black curve corresponds to the singlino-dominated 
flat part of the S- model in Figure El with M 2 < 0. Fixed values of h s = 0.75 and s = 250 GeV are 
used. The acceptable relic density is VL x oh 2 ~ 0.1 . 

• When m x o ~ Mz/2, there is a dip in the curve of relic density versus m x o since there 
is a peak in the annihilation cross section enhanced by the Z-poIe. 

• In the S'-model, a positive value of M 2 gives two acceptable CDM density regions 
(around the Z pole and the enhanced coupling region such as m x o ~ 80 GeV in the 
tan/3 = 1.03 case) while a negative M 2 leads to Q x oh 2 > 1 except for the singlino- 
dominated region. 

• The MSSM can also give an acceptable CDM density, but the model is excluded for 
tan/3 ~ 1 by the LEP Higgs bound, owing to its small Higgs mass. The S'-model 
Higgs contains a singlet component and the light Higgs can be compatible with the 
LEP constraint [3jJ. 

• The small isolated region for relatively large m x o corresponds to M 2 < with the 
lightest neutralino being singlino-dominated in the S- model 12 . As observed in Figure 
01 a sudden sizable deviation of |-/V 13 | 2 from \Ni^\ 2 in a singlino dominated region 
(M 2 < —150 GeV for tan (3 = 1.03) provides a sudden drop of Q x oh 2 from the rest of 
the M 2 < curve. 

Figure 03 presents ranges of neutralino relic density in regions of the M 2 -s plane in the 
S-model. We choose h s = 0.75 with (a) tan/3 = 1.03, (b) tan/3 = 2.0, and (c) tan/3 = 2.5. 
Also shown is the region excluded by the LEP 2 chargino mass bound with m x + < 104 GeV 
[36| The parameter points at which the light chargino is the LSP were omitted. The tree 
level mass for the lighter chargino {m x + < m x + ) is evaluated with the chargino mass matrix 

/ M 2 V2M w sm(3\ 



12 The case in which the lightest neutralino is mostly singlino-like (e.g., tan/3 = 1.03, m x o ~ 80 GeV for 
both M-2 < and M 2 > 0), is qualitatively similar to those studied in Ref. Q. 
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FIG. 5: Regions of the S'-model neutralino relic density in the M2-S plane (with s scanned only 
above 100 GeV) for 0.09 < tt x oh 2 < 0.15 (filled square; 3a allowed range), n x oh 2 < 0.09 (open 
circle), and 0.15 < £l x ah 2 < 1.0 (cross). Three representative values of tan/3 are chosen: (a) 
tan/3 = 1.03, (b) tan/3 = 2.0, and (c) tan/3 = 2.5, and a fixed value of h s = 0.75 is used. The 
shaded region of the parameter space (bounded by solid curves) is excluded by the LEP 2 chargino 
mass limit (m + ^$ 104 GeV). There exist sizable regions (filled square) in the parameter space 
consistent with the relic density constraint outside of the chargino mass exclusion boundary. 

We present the S'-model neutralino relic density in 3 regions: 0.09 < Q x oh 2 < 0.15 (filled 
square), Q x oh 2 < 0.09 (open circle), and 0.15 < Q x oh 2 < 1.0 (cross). The 3a range 13 of the 
CDM relic density of Eq. (JJJ is 0.09 < Q x oh 2 < 0.15; however if there are other sources of 
dark matter in addition to the lightest neutralino, the range Q x oh 2 < 0.09 would be relevant. 
Due to the finite grid, the filled square should be understood to be on the boundary of the 
open circle and the cross. 

There are three separate regions that have an acceptable CDM density (filled square): 
(PI) Near the Z pole, (P2) enhanced coupling region (tan/? « 1 case), and (P3) isolated 
singlino region. In general, the Zx°X° coupling is enhanced by the singlino component. 
There appears a sudden peak of the enhancement (P2) when tan /3 ~ 1 as we can see from 
Figure El (a). In the enhanced coupling region (P2), even when the x° mass is significantly 
distant from the Z pole, an acceptable relic density can be obtained. The isolated singlino 
region (P3) is singlino-dominated and happens for M 2 < 0. 

For tan/3 = 1.03, there is a small region in the M 2 -s plane with M 2 > that satisfies the 
relic density and LEP chargino mass constraints. The solution in this region is due to the 
enhanced Zx°X° coupling. For tan/3 = 2.0, there is a large acceptable region with M 2 < 0. 
With M 2 > most of the parameter regions that give the relic density are excluded by the 
LEP 2 chargino search. For tan/3 = 2.5, there is a large region with M 2 > that reproduces 
the observable relic density and is consistent with the chargino mass limits. 

We numerically checked, with suitable parameter values including different h s , that the 
S'-model can reproduce the observed relic density for most of the theoretically allowed m x o 
range without violating the LEP constraints. However, for a relatively large neutralino 
mass, i.e., m x o w 80 to 100 GeV, it becomes hard to satisfy both the relic density and the 



Since we are using tree level masses and couplings, we allow rather conservative 3<r range for the allowed 
CDM relic density. 
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LEP constraints since those m x o values occur only with tan (5 close to 1 and h s close to its 
maximum (0.75), but then the chargino constraint becomes severe, as can be seen in Figure 
El (a). 

IV. CONCLUSIONS 

We have studied the properties of th e lig htest neutralino in a supersymmetric U(l)' model 
with a secluded U(l)' breaking sector [25j. In this model, tan/3 ~ 1 is required for natural 
electroweak symmetry breaking. In general, the model allows a lightest neutralino mass up 
to about 300 GeV (depending on the charge assignments and input mass limits) with the 

2 

gaugino unification assumption M\> = Mi = |%M 2 , but only up to m x o < 100 GeV in 

the limit with My 3> Mi = |^jM 2 . We quantitatively studied the S'-model in the limit 

where Mi , Sj are much larger than the electroweak scale. In this limit the Zx°X° coupling 
is enhanced compared to the MSSM due to the singlino component. This allows the x° P arr 
annihilation rate via the Z resonance channel alone to reproduce the observed relic density 
of the cold dark matter. In addition, the doublet-singlet mixed nature of the Higgs bosons 
allows the lightest Higgs bosons to have a small mass without violating the LEP constraint 
on SM-like Higgs mass. The S'-model explains the relic density over a considerable fraction 
of the parameter space. 
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